Role of the exchange and correlation potential into calculating the x-ray absorption 
spectra of half-metallic alloys: the case of Mn and Cu K-edge XANES in Cu 2 MnM 

(M = Al, Sn, In) Heusler alloys. 
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This work reports a theoretical study of the x-ray absorption near-edge structure spectra at 
both the Cu and the Mn K-edge in several Cu2MnM (M= Al, Sn and In) Heusler alloys. Our 
results show that ab-initio single-channel multiple-scattering calculations are able of reproducing 
the experimental spectra. Moreover, an extensive discussion is presented concerning the role of the 
final state potential needed to reproduce the experimental data of these half-metallic alloys. In 
particular, the effects of the cluster-size and of the exchange and correlation potential needed in 
reproducing all the experimental XANES features are discussed. 

PACS numbers: 78.70Dm, 61.10Ht 
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I. INTRODUCTION 

X-ray absorption spectroscopy (XAS) constitutes 
nowadays an outstanding structural tool. The extended 
x-ray absorption fine structure (EXAFS) part of the spec- 
trum is commonly used to determine the local environ- 
ment around a selected atomic species in a great vari- 
ety of systems^. Moreover, the near-edge part of the 
absorption spectrum (XANES) becomes an incompara- 
ble stereochemical probe because its high sensitivity to 
the bonding geometry. This capability is due to the low 
kinetic energy of the photoelectron that favors the con- 
tribution of multiple scattering processes. By this rea- 
son, a great effort has been devoted in the last decade to 
obtain structural determinations from XANES, includ- 
ing bond-angles information. However, the interpreta- 
tion and the ab-initio computation of XANES spectra 
are not so straightforward as for EXAFS. 

XANES computation requires sophisticated simula- 
tion tools 2 -^^! 7 -^. Most of these ab-initio codes are 
built within the one-electron multiple-scattering (MS) 
framework 9 - and by using the so called muffin-tin approx- 
imation. More recently, several works appear reporting 
XANES computations performed within a non-muffin-tin 
approac h 7 : 10 ' 11 and the promising implementation of the 
multi-channel MS theor y 12 : 13 . Despite the above pro- 
gresses, the construction of the scattering potential still 
remains an open problem. Typically, the standard ab 
initio XANES calculations are performed in the frame- 
work of the muffin-tin approximation with the addition 
of an exchange-correlation potential (ECP). In the case 
of complex ECP potentials, the imaginary part also ac- 
counts for the damping of the excited photoelectron. 
At present, the energy dependent Hedin-Lundqvist (HL) 
complex potentia l 14 : 15 is the most widely used into ab 
initio XANES calculations^ 6 -. However, it has been re- 
ported that in several cases the energy dependent Dirac- 
Hara (DH) exchange potential shows the best agreement 
to the experimental dat a 17 : 18 ' 19 : 20 , although this assign- 



ment is not free of controvers y 22 : 23 : 24 : 25 . As a result, the 
choice of the ECP among the usual scattering potentials 
set, ranging from real X Q to complex Hedin-Lundqvist, 
becomes one the most important steps into obtaining a 
good reproduction of the experimental spectra. The ex- 
perience accumulated through the calculation of the ab- 
sorption cross section for different systems suggests that 
the HL-ECP offers a good performance in the case of 
metals, while for ionic and covalent systems the use of 
DH seems to be more adequat o 17 : 18 ' 19 : 20 ' 21 . However, 
the variety of reported results indicates that this choice 
is material specific and, unfortunately, prevents one of 
a-priori fixing the ECP for each particular case. In this 
work we have faced this problem in connection to the 
case of half- metallic alloys. These systems present a pe- 
culiar behavior associated to the metal and insulator or 
semiconductor character of the spin subsystems. 

The Cu 2 MnM (M= Al, Sn and In) Heusler alloys are 
interesting magnetic systems because they possess lo- 
calized magnetic moments although they are all metal- 
lic. The magnetic, electrical and structural properties 
in these systems are known to be strongly dependent 
on both the conduction electron concentration and the 
chemical orde r 26 : 27 ' 28 : 29 : 30 . Despite the importance of 
getting a deeper insight in the electronic structure and 
the chemical bonding mechanism in these compounds, 
few amount of research has been performed regarding the 
magnetic polarization of the conduction electrons and its 
role into the magnetic ordering. Recently, Uemura et al. 
have reported the results of a combined XAS and x-ray 
magnetic circular dichroism (XMCD) study performed 
at the K-absorption edge of both Mn and Cu in these 
Cu 2 MnM Heusler alloys^. Theoretical calculations have 
been made on the basis of photoelectron full multiple 
scattering (FMS) theory^ 2 -, and by using the Madclung 
potential taken from the band calculation for Cu^MnAl 
performed by Ishida et a/4 -. According to Uemura et 
al. the spin-orbit interaction and the exchange scatter- 
ing are the main factors for the calculation. Thus, even 
when they used a non-local exchange potential that can 
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directly take the spin-polarization effects into account, 
correlation-potential was not included. In addition, only 
51 atoms, 5 A, were included to built up the cluster used 
for the computations. The fact that computations were 
performed for such a small cluster and without correla- 
tion potential poses serious concerns regarding the relia- 
bility of the reached conclusions. 

In this work, we present detailed ab initio computa- 
tions of the Mn and Cu K-edge XANES spectra in the 
Cu 2 MnM (M= Al, Sn and In) Heusler alloys. These 
calculations of the absorption spectra have been made 
within the multiple-scattering framework. Special atten- 
tion has been paid to establish the improvements ob- 
tained by using different treatments of the exchange- 
correlation part in the final state potential, intending to 
obtain the best agreement both in energy and intensity, 
between the experimental and theoretical XANES spec- 
tra. 



II. COMPUTATIONAL METHODS 

The computation of the XANES spectra was carried 
out using the multiple-scattering code CONTINUUM^ 
based on the one-electron full-multiple-scattering 
theoryiLLS, A complete discussion of the procedure can 
be found in Ref^ and Refi^ 3 .. Computations were made 
in parallel mode by using the MPI library 3 ^. 

The potential for the different atomic clusters was 
approximated by a set of spherically averaged muffin- 
tin (MT) potentials built by following the standard 
Mattheis' prescription 35 . The muffin-tin radii were de- 
termined following the Norman's criterion and by impos- 
ing an overlapping factor ranging from 1% to 10% 36 . We 
have verified that the factor chosen within this range does 
not affect the spectral shape and only slight differences 
are found in the first 5 eV of the computed spectra. Af- 
ter convolution of the spectra to account for both the 
core-hole lifetime and the experimental resolution these 
differences are not significant. 

The Coulomb part of each atomic potential was gener- 
ated using charge densities from atomic code of non-local 
selfconsistent Dirac-Fock code 5 -^. The atomic orbitals 
were chosen to be neutral for the ground state potential. 
Two different approximations were tested for the excited 
state potential: i) the same potential as for the ground 
state was used; ii) relaxed Z+l approximation 38 . During 
the present calculations we have found that the screened 
and relaxed Z+l option leads to the best performance 
into simulating the experimental absorption spectra at 
both the Mn and Cu K-edge. Finally, we have tested 
three different choices for the exchange and correlation 
part of the final state potential: X Q , the energy depen- 
dent Hedin-Lundqvist (HL) complex potential and the 
energy-dependent Dirac-Hara (DH) exchange potential. 

The computed spectra have been compared to the ex- 
perimental XANES data reported Uemura et al£^. In 
all the cases, the calculated theoretical spectra have been 



FIG. 1: The stoichiometric Cu 2 MnM unit cell: Cu (black 
circles), Mn (crossed circles) and M (open circles). 



further convoluted with a Lorentzian shape function T= 
1.5 eV to account for both the core- hole lifetime^ and 
the experimental resolution. 



III. RESULTS AND DISCUSSION 



Ab-initio XANES calculations have been performed at 
both the Mn and Cu K-edge in the case of Cu2MnAl, 
Cu2MnIn and Cu2MnSn. These Heusler alloys show the 
L2i-type (space group FmSm) ordered structure illus- 
trated in Fig. [T] The unit cell contains eight of the small 
body-centred cubic cells that can be regarded as four 
interpenetrating fee sublattices. There are four kinds 
of atomic positions: Cu atoms occupy (1/4, 1/4, 1/4) 
and (3/4, 3/4, 3/4) positions, Mn atoms are at (1/2, 
1/2, 1/2), while the M (Al, In, Sn) atoms occupy (0, 0, 
0) position. The lattice parameter is a = 5.957 A for 
Cu 2 MnAl- 40 ; a = 6.166 A for Cu 2 MnSr4i; and a = 6.217 
A for Cu 2 MnIn 4 ^. 

First step of the calculations was deserved to deter- 
mine the size of the cluster needed to reproduce all the 
spectral structures present in the experimental XANES 
spectra. To this end we have performed the computa- 
tion of the Mn K-edge XANES in the case of Cu 2 MnAl 
by increasing progressively the number of atoms in the 
built-up cluster. In this way we have computed the Mn 
K-edge absorption for clusters containing: 59 atoms, i.e. 
including contributions from neighbouring atoms located 
within the first 5.2 A around photoabsorbing Mn; 65 (6 
A), 113 (6.7 A) and 137 (7.3 A) atoms. The compu- 
tations, shown in Fig. [2J have been performed by using 
the real part of the HL ECP and with lm ax = 4. The 
computed spectra show the dependence of both the en- 
ergy position and the relative intensity of the spectral 
features as the size of the cluster increases. In the case of 
the smallest clusters, 59 and 65 atoms, several spectral 
features are not reproduced by the calculation. In partic- 
ular, the B 2 peak ( E ~ 26 eV), located at the high energy 
side of the main Bi peak is not reproduced. Moreover, 
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FIG. 2: (color online) Comparison between the experimental 
Mn K-edge of Cu2MnAl (•) and the MS computations per- 
formed for different cluster sizes covering up to 7.3 A around 
the absorbing Mn: 59 (black, solid line), 65 (blue, dotted 
line), 113 (red, dashed line) and 137 atoms (green, o). 

for these clusters the computation yields a broad reso- 
nance at ~ 50 eV, whilst the experimental spectrum ex- 
hibits two well defined peaks Di and D2. By contrast, the 
computations performed for bigger clusters, 113 and 137 
atoms, account for the previously missed structures. In 
addition, no differences are found within the first 60 eV 
of the absorption spectrum between computations per- 
formed for clusters including contributions from atoms 
located within the first 6.7 A and 7.3 A around the pho- 
toabsorber. This result indicates that the addition of fur- 
ther coordination shells does not contribute significantly 
to the XANES spectrum. Accordingly, all the calcula- 
tions reported henceforth have been obtained by using 
169 atoms to generate the scattering potential but only 
the scattering contributions of the first 137 atoms (~ 6.7 
A) are computed. 



Next, we have investigated the maximum angular mo- 
mentum quantum number, l ma x, needed to account for 
the experimental spectrum in the first 100 eV of the ab- 
sorption spectrum. The choice of l max = 3 and 4 does 
not affect the result of the computations in the first 30 
eV of the spectrum above the edge. However, for higher 
energies both the energy position and intensity of the 
computed features clearly differ. By increasing l ma x from 
4 to 5, the computed spectra do not change in the sim- 
ulated spectra in the energy interval of interest, up to ~ 
80 eV above the edge. Therefore, we have fixed lmax = 
4 for all the calculations. 

Based on the above results we have fixed both the clus- 
ter size and lmax- Then, we have investigated the role of 



tjJj^ 2different exchange and correlation potentials (ECP) 
into correctly reproduce the experimental spectra. We 
have computed the Mn K-edge XANES of Cu 2 MnAl by 
using the Xa and the energy dependent Hcdin-Lundqvist 
(HL) and Dirac-Hara (DH) ECP potentials^. The HL 
is a complex potential in which its imaginary accounts 
for the inelastic losses of the photoelectron. Therefore, 
we have also made computations by using only the real 
part of the HL ECP (hereafter, real HL) and, in addition, 
we have built a "complex" DH one by adding to it the 
imaginary part of the HL one. The results of the com- 
putations are shown in Fig. [31 For the sake of clarity, 
the comparison between the theoretical spectra and the 
experimental one has been divided in two panels, one for 
the real potentials (a) and the second for the complex 
ones (b). In the case of the real potentials, the inten- 
sity of the spectral features is not in agreement with the 
experimentally observed, as expected because the pho- 
toelectron damping is not taken into account. However, 
its inspection is rather useful as it allows one to deter- 
mine that the computations reproduce all the spectral 
features shown by the experimental spectrum. To this 
respect, the agreement between the computations and 
the experimental signal is noticeable as the calculations 
reproduce the complex profile at the threshold (Ai and 
A2 features), the main peaks (Bi, Di, D2 and F) and the 
two tiny structures, B2 and C. The performance of the 
different potentials is similar although in the case of both 
Xa and the real HL the calculated absorption maxima 
falling short of the observed ones. The addition of the 
imaginary part of the HL potential improves the calcu- 
lations. Good agreement between the computations and 
the experimental spectrum is obtained for both complex 
Hedin-Lundqvist ECP and by adding to the Dirac-Hara 
(DH) exchange the imaginary part of the HL (hereafter 
complex Dirac-Hara). As shown in Fig.[3jD) all the exper- 
imental features are properly accounted by the compu- 
tations. Complex HL computation yields a better repro- 
duction of the intensity ratio among the different spectral 
features, specially regarding the main peak Bi and both 
Ai and B2 features. By contrast complex DH shows a 
better reproduction of the relative energy position among 
the peaks. The experimental energy separation between 
the main absorption peak Bi and the negative E deep is 
AE = 48 eV, while computation yields 47.6 eV and 46.2 
for complex DH and HL, respectively. 



We have performed the same class of computations for 
both Cu2MnIn and Cu2MnSn systems. In all cases we 
have used lmax = 4 and the potential was created by 
using a 169 atoms cluster to calculate the scattering con- 
tributions of the first 137 atoms. In this way, we ac- 
count for the contributions from atoms located within 
the first 6.9 A and 6.95 A around Mn for Cu 2 MnSn and 
Cu2MnIn, respectively. As shown in Fig. Q] good agree- 
ment between the ab-initio calculations and the experi- 
mental data is also found for these compounds. Similarly 



4 




-20 20 40 60 80 100 

E - E (eV) 



FIG. 3: (color online) Comparison between the experimen- 
tal Mn K-edge of Cu2MnAl (o) and the MS computations 
performed for different ECP potentials: a) Xa (green dashed 
line), real HL (red, solid line) and real DH (blue, dotted line); 
b) complex HL (red, solid line) and complex DH (blue, dotted 
line) (see text for details). 



to the CU2M11AI case, the use of the complex DH poten- 
tial shows a better performance that the complex HL 
one into reproducing the energy difference, AE, among 
the spectral features. Experimentally, the energy differ- 
ence between the main peak (B) and the first negative 
deep (E) is 26 eV for both Cu 2 MnIn and Cu 2 MnSn com- 
pounds. HL computations yield AE = 23.3 and 22.5 eV 
for the In and Sn compounds, respectively. These results 
are improved by using the complex DH ECP yielding <5E 
= 25.9 eV and 24.8 eV for Cu 2 MnIn and Cu 2 MnSn, re- 
spectively. As observed for Cu 2 MnAl, the intensity ratio 
of the different spectral features in the first 25 eV's of 
the spectrum is better reproduced by using HL than DH 
ECP. Indeed, the weak structure (C) at the high energy 
side of the main peak (B) is correctly reproduced by the 
HL computation. By contrast, the absorption maxima 



calculated by using the HL ECP fall short the experimen- 
t'&l 3 ones, specially regarding the energy region between 30 
to 70 eV above the edge, (D to G structures, Fig.S]). 

We have extended our study to the Cu K-edge in the 
same materials by following identical procedure as for the 
Mn K-edge. Also in this case, we have verified the need 
of using Imax = 4 as well as big clusters to reproduce 
the Cu K-edge experimental spectra. Best agreement 
between the experimental data and the calculations have 
been obtained by using 169 atoms to generate the poten- 
tial and taken into account the scattering contributions of 
the first 137 atoms. Therefore, contributions from atoms 
located within the first 6.7 A, 6.9 A and 6.95 A around 
Cu are taken into account for Cu 2 MnAl, Cu 2 MnSn and 
Cu 2 MnIn, respectively. 

The experimental Cu K-edge XANES spectra are sim- 
ilar in the three cases. There are two contributions at 
the threshold (Ai and A 2 ) followed by a prominent peak 
(B) at ~ 15 eV in the three compounds. In the case 
of Cu 2 MnAl, there is a single broad resonance (D) at 
about 27 eV above the edge. By contrast, the spectral 
profile in this region is different for both Cu 2 MnSn and 
Cu 2 MnIn compounds. The Cu K-edge XANES of these 
compounds shows a small structure (C) on the low en- 
ergy side of the broad one (D). The emerging feature is 
located at ~ 23 eV while the broad D one is shifted to 
higher energies (~ 32 eV) with respect to the Cu 2 MnAl 
case. The theoretical calculations are able of reproduc- 
ing this trend. As shown in Fig. [5] the computed spectra 
for both Cu 2 MnSn and Cu 2 MnIn clearly differ from that 
obtained for Cu 2 MnAl in the energy region correspond- 
ing to C and D features, i.e. from 20 to 35 eV above 
the edge. The best reproduction of the Cu K-edge ex- 
perimental spectra is obtained by using the complex DH 
potential. Contrary to the Mn K-edge, the improvement 
of the complex DH ECP over the HL one concerns both 
the intensity ratio and the relative energy separation be- 
tween the different spectral features. Indeed, by using 
the complex HL ECP the relative intensity of the main 
contributions to the absorption main peak is not satisfac- 
tory and the contraction of the energy scale still persists. 



The fact that the agreement between the theoretical 
calculations and the experiment is better when using the 
complex DH instead of HL potential is noticeable. In- 
deed, it is commonly assumed that the HL potential, 
with its energy dependent exchange and its imaginary 
part, is able to reproduce rather accurately the experi- 
mental spectra in the case of metals. By contrast, sev- 
eral works have reported the improved performance of 
the DH potential into accounting for the experimental 
spectra of ionic systems, insulators and transition metal 
oxide o 17 ' 18 ! 19 ' 24 . This result suggests that the treat- 
ment of the Coulomb hole correlation by using a single 
plasmon-pole within the GW approximation 4 ^ is rather 
crude in the high energy region. HL potential is sepa- 
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FIG. 4: (color online) Comparison between the MS compu- 
tations performed by using both complex HL (red, solid line) 
and complex DH (blue, dotted line) and the experimental Mn 
K-edge (o) of: a) Cu 2 MnIn and b) Cu 2 MnSn. 
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rated into three parts, 

V HL = V° H + V S ex + V ch (1) 

where Vj^ H is the bare exchange potential, Dirac-Hara, 
V sex is the screening of exchange potential and V c h is 
the Coulomb correlation hol o 14 ' 15 . Vj^ H is negative and 
at high energy asymptotically behaves like, fc~ 2 . Ac- 
cording to several authors this energy dependence is cru- 
cial in the high-energy part of the absorption spectrum, 
EXAF S 44 ' 45 . V sex is positive and roughly constant and 
then rapidly decreases like k , while V c h, which de- 
scribes the repulsion of electrons of the same spin, is 
negative and its behavior is fc _1 at high energy^. As 
discussed in Ref4^ Coulomb correlation hole will lower 
the energy of the states and, as a consequence, the po- 
sition of the peaks in continuum state is lowered. The 
Coulomb correlation hole effect dominates in the higher 
energy region, in agreement to our results. Indeed, the 



FIG. 5: (color online) Comparison between the 
experimental(o) Cu K-edge of (a) Cu2MnAl, (b) Cu2MnIn 
and (c) Cu2MnSn, and the result of the MS computations 
performed for both complex HL (red, solid line) and complex 
DH (blue, dotted line) ECP potentials. 



comparison of the computations performed by using DH 
and HL ECP shows that in the latter case all the peaks 
in the high energy region are shifted to a lower energy 
position. This result reflects the inadequacy of using a 
single plasmon-pole within the GW approximation 4 ^ to 
account for the Coulomb hole correlation in the high en- 
ergy region. 
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IV. SUMMARY AND CONCLUSIONS 

We have presented the detailed ab-initio computation 
of the Mn and Cu K-edge XANES spectra in the case 
of Cv^MnM (M= Al, Sn and In) compounds performed 
within the multiple-scattering framework. 

The comparison between the experimental data at the 
two edges and the theoretical calculations has demon- 
strated the need of using large clusters, ~ 7A around 
central photoabsorbing atom, either Mn or Cu, to repro- 
duce the experimental spectra of these Heusler alloys. 

During the multiple-scattering ab initio calculations we 
have tested different exchange and correlation possibili- 
ties to construct the final state potentials. Contrary to 
the current expectation for metallic systems of obtain- 
ing the best reproduction of the experimental XANES 
spectra by using the Hedin-Lundqvist potential, the best 
agreement is obtained, specially in the Cu K-edge case, 
by using the Dirac-Hara ECP. 

We concluded that Hedin-Lundqvist potential overes- 
timates the effect of the Coulomb correlation hole in the 
high energy region resulting in an unphysical contraction 
of the computed spectra by lowering in energy the calcu- 
lated peak positions. This problem is solved by adding 



to the real energy-dependent Dirac-Hara potential the 
imaginary part of the Hedin-Lundqvist ECP. In this way, 
the inelastic losses of the photoelectron are accounted 
and both the relative energy position and the intensity 
ratio among the different spectral features are properly 
accounted by the computations. 
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